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Abstract 

A recent model of hadron states is extended to include meson decays. We find 
that the overall success of the model is quite good. Possible improvements to 
the model are suggested. 



I. INTRODUCTION AND MOTIVATION 



Recently, a model for the description of hadronic states has been developed and applied 
to both meson and baryon spectra JI|fJ] . The success of the model in describing the spectra in 
these sectors has led to attempts to take the model further by applying it to multiquark states 
0. This success has also led to the question of whether the model describes other aspects of 
baryon and meson spectroscopy adequately. In particular, it is useful to know whether the 
model can help shed light on some of the many outstanding puzzles in hadron spectroscopy. 
Before such a question can be addressed, the model should be able to reproduce those 
portions of hadron phenomenology that are better understood. In this article, we present 
the results of application of the wave functions obtained from this model to the strong 
decays of light mesons. We find that the model is as successful as other models, for the 
decays examined. 

The rest of this article is organized as follows. In the next section we summarize the 
salient points of the hadronic model, as well as the model we use to describe the decays of 
hadrons. In section 3 we present our results with some discussion, while our conclusions are 
presented in section 4. 

II. THE MODEL 

Two ingredients are essential for a good description of meson decays. One of these is 
a realistic description of the mesons themselves. The second is a valid description of the 
transition operator. Both of these depend very much on our educated guesses about non- 
perturbative QCD dynamics. In the case of the former, there are many models proposed. 
We confine ourselves to the work of [|TJ. For the meson decays, there are also a number of 
models proposed, but we will limit ourselves to the pair creation model first proposed by 
Micu {|, and popularized by the Orsay group ||. 

A. Spectrum 

One possibility for the analysis of the meson spectrum is the use of a 'good' quark- 
antiquark potential. This is the course we pursue. Although the one-gluon exchange mech- 
anism suggests possible functional forms for some terms of the potential, the fact that we 
still do not know much about the mechanism of confinement means that we must rely on 
phenomenology to determine both the form and the parameters of the potential. Several 
potentials have been proposed in the past in the literature. We believe that in models of the 
kind that we propose, a very crucial test is a unified description of the meson and baryon 
spectra. This has been undertaken by one of us |],[2|], and six different potentials have been 
found which fulfill this criterion. These potentials have been applied, not only to the meson 
and baryon sectors, but also to tetra-quark states || , as well as to a description of the K N 
interaction ||. The results in all of these cases have been encouraging. 

In this paper, we select two of these potentials for this study. The general form of each 
potential has been more or less imposed by some basic QCD constraints, but the parameters 
have been determined by a fit to a well-chosen sample of meson and baryon states. Both 
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potentials are used with a non-relativistic kinetic energy term, and they take the general 
form 




- + Ar p + A + 



2ttk' exp(— r 2 /rl) 




(1) 
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Smiirij 7r 3 / 2 7"o 



One peculiarity of these potentials is that the range r of the hyperfme term is mass 
dependent through the relation 



In the limit of a vanishing range, this term gives the usual Fermi-Breit prescription 8{f). 

These potentials contain neither spin-orbit nor tensor contributions so that 3 Lj states 
have the same wave functions. In the same spirit, 3 Lj <-> 3 L'j and 3 Lj ^ x Lj mixings are 
absent. This approximation can have some effect on the results, but we are mainly interested 
in the gross features of the decays, so that the potentials are quite satisfactory from this 
point of view. Furthermore, there is no isospin dependence (this is a property shared by 
most potential models: in M , isospin dependence is introduced by means of mixing through 
annihilation). Thus, the p and u, as well as the aj and fj resonances are degenerate, both 
in energy and in wave function. More dramatic is the case of the pseudoscalar mesons, as 
this leads to a degeneracy of the 7r and the non-strange part of the 77. Isospin dependence 
appears only through the flavor part of the wave function, although there probably should 
be a difference in the spatial part of the wave function as well. In any case, we consider two 
potentials that yield good overall results for the meson spectrum. 

The first potential, denoted AL1, has the usual Coulomb+linear form for the central 
part, and its parameters are 



The second choice of potential, denoted API, has a confining term suggested by the 
Regge trajectory behavior of orbital states in a non-relativistic treatment. The parameters 
are 



Both potentials reproduce spectra of comparable quality and are simple enough to be 
handled without any difficulty. In particular, it is quite easy to solve the differential equation 
resulting from the Schrodinger equation. However, the radial part R n i. s (r) of the meson wave 
function is calculated numerically on a grid. Such a form is not easily used for studying 
more complicated problems: a continuous expression is usually more convenient. 

To solve this problem, we approximate the regularized part of the exact radial wave 
function by a linear combination of gaussian functions, 




(2) 



m u = m d = 0.315GeV; m s = 0.577GeV; m c = 1.836GeV; m b = 5.227GeV 
re = 0.5069; re' = 1.8609; A = 0.1653GeV 2 ; p = 1; A = -0.8321GeV 
B = 0.2204; A = 1.6553GeV B_1 . 



m u = md = 0.277GeV; m s = 0.553GeV; m c = 1.819GeV; m b = 5.206GeV 
re = 0.4242; re' = 1.8025; A = 0.3898GeV 5/3 ; p = 2/3; A = -1.1313GeV 
B = 0.3263; A = 1.5296GeV B_1 . 
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TABLE I. Masses of some mesons used in the paper obtained with potential ALL The exper- 
imental values are given in the second column; for the aj and fj resonances we show the energy 
centroid. The third column reports the exact values obtained from the Schrodinger equation by 
solving the resulting differential equation. The last three columns deal with the approximation 
based on gaussian functions with, respectively, N = 1, 2, 3 terms in the expansion. 



N 

Rnis(r) = £ c t r l exp(-a,r 2 /2). (3) 

For a given number N of gaussian functions, the parameters Cj and are determined by 
a variational procedure on the energy of the considered state. = 1 is a rather rough 
approximation, but A^ = 2 or A^ = 3 greatly improves the results. This is illustrated in 
figure [I], where the exact wave function, as well as the A^ = 1,2,3 approximations, are 
plotted for the pion, for a large range of r values. In this figure, the curves on the left are 
for the potential AL1, while those on the right are for API. We see that except for very 
short distances, the approximation with 3 gaussian terms gives essentially the exact wave 
function. This expression, eqn. ([!]), is then used in our model for meson decays. 

We show in table | a few of the masses obtained in this way, for the potential All. 
Experimental values are shown in the second column, while the exact predictions are given 
in the third. The masses obtained using the expansion of the wave function in gaussian 
functions are shown in the next three columns for AT = 1, 2, 3, respectively. We see that the 
approximation with A^ = 3 gives essentially the same results as the exact case, confirming 
the trend illustrated in figure 0. 

B. Decays 

For the OZI || allowed decays of a meson, we assume that quark pair creation, somewhere 
within the volume of the parent hadron, provides the dominant contribution. In fact, for the 
operator we choose, this creation can take place anywhere in space, but the wave function 
overlaps ensure that the dominant contribution to the amplitude comes from creation within 
the 'confines' of the parent hadron. 

The full mechanism for strong decays is, of course, as yet unknown, but the success of 
pair creation models hints at what may be the most important aspects of the underlying 
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FIG. 1. The pion wave function that results from two choices of potential, with different ex- 
pansion approximations. The curves on the left are from potential AL1, while those on the right 
are from API. 




dynamics. We assume that the operator responsible for the decay is 
T = -37 J2 J dpidpjdipi+p^CijFij 

x £ < l,m;l, -m\0,0> X™yi m (ti-Pj)b\(Pi)d f M)- (4) 



Here, and are the color and flavor wave functions of the created pair, both assumed 
to be singlet, Xij is the spin triplet wave function of the pair, and 3MPi — Pj) is the vector 
harmonic indicating that the pair is in a relative P-wave. b\ creates a quark, and d\ creates 
an antiquark. 

For the transition A — > BC, we are interested in evaluating the transition amplitude M, 
given by 

M=<BC\T\A>. (5) 

The details of this calculation are given elsewhere HfjlJ. Once M is evaluated, we calculate 
the decay rate as 

T(A^BC)= ^ \M A ^ BC (J bc ,£,k )\ 2 ^(ABC), (6) 

Jbc/ 

where $ is the phase space for the decay. We have investigated three possible choices for 
this phase space. The first choice is that of fully relativistic phase space, for which 

= 27r E b (ko)E c (k )k 
m a 



with Ei(ko) = y&o + mj, and k is calculated in the rest frame of the parent meson. The 
second choice is fully non-relativistic phase space, with 
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$(ABC) = 2n^^. (8) 



Our third choice is the 'weak-binding' phase space suggested by Kokoski and Isgur [yj, for 
which 

$(ABC) = 27;^^, (9) 
m a 

where the m's are effective meson masses (in the model) without hyperfine interaction. For 
the rhi, we use the values suggested in |TT|. 



III. RESULTS AND DISCUSSION 

We have performed a number of fits to the experimental data; we show the results of 
these fits in table ||. In addition, using the parameters obtained from these fits, we have 
calculated the amplitudes for a number of other decays. These are presented in tables [TTT|, [TV] 
and [V|. In table U, we consider only decays in which the daughter mesons are narrow. For 
this purpose, we treat the us as a narrow state. In tables [TTT] to |V|, we extend the calculation 
to include p's and K*'s as daughter mesons. 

The numbers in column two (Fit I) of these tables result from a fit to the amplitude for 
p — ► 7T7T, using the usual form of the 3 Po operator. The numbers in column three (Fit II) 
result from using the same operator, but fitting to the 22 decays shown in table |I[ Column 
four (Fit III) corresponds to a similar global fit, but using a modified transition operator. 
In this case, the operator has been modified by making the change 

lY,X m yr m (P) <l,m;l,-m|0,0> 

m 

- e"^ (71 + I2V 2 ) E X^r™ (p) < 1, m; 1, -m\0, >, (10) 

m 

with A = 1.37, 71 = 0.14 and 72 = 1.4. For the results of column 2 and 3, we find 7 = 0.41 
and 7 = 0.37, respectively. These three fits all result from the AL1 wavefunctions, expanded 
to N = 1. 

The effects of the size of the expansion basis of the wave functions are illustrated in Fits IV 
and V; the former is a one parameter global fit (7 = 0.31) using the AL1 wave functions with 
three gaussians, while the latter uses the same wave functions with the modified transition 
operator (A = 0.0, 71 = 0.14, 72 = 0.94). Fits VI and VII show the effects of a change 
in potential. Both of these arise from the API wave functions in conjunction with the 
modified transition operator: Fit VI uses wave functions with a single gaussian (A = 1.1, 
71 = 0.13, 72 = 1.3), while the wave functions of Fit VII are expanded in a basis of three 
gaussians (A = 0.0, 71 = 0.10, 72 = 1.0). For Fits III, V, VI and VII, the theoretical 
errors are evaluated using the covariance matrix resulting from the fit. All of the amplitudes 
are calculated using the weak-binding prescription for phase space. As can be seen from 
tables 0, there are few instances in which a change of the choice of potential or wave 
function expansion, accompanied by apparently radical changes in fit parameters, induces 
a big change in the amplitude. In particular, the modification of the operator induces an 
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TABLE II. Results from various fits. Fit I corresponds to a fit to the single decay p — > tttt, 
using the simple form of the 3 Po operator. Fit II corresponds to the same choice of operator, 
but to a global fit to all of the decays in this table. Fit III also results from a global fit, using 
the modified form of the 3 Po operator. Fits I, II and III result from single-component, AL1 wave 
functions. Fits IV to VII are all global fits to all the decays shown in the table. Fit IV results from 
using the three-component, AL1 wave functions, with the usual form of the 3 Po operator. Fit V 
uses the same wave functions and the modified 3 Pq operator. Fit VI and VII also use this form 
of the operator, but Fit VI uses the single component form of the API wave functions, while Fit 
VII uses the three-component form of the API wave functions. All widths are calculated using the 



weak-binding phase space prescription. The experimental results are from [12 



decay 












r x /2 (Mcv 1 ^) 












Fit I 


Fit II 


Fit III 


Fit IV 


Fit V 


Fit VI 


Fit VII 


experiment 


p — y 7T7T 




12.3 


11.1 


11.9 ±0.3 


10.4 


12.1 ± 0.3 


12.0 ± 0.3 


12.2 ± 0.3 


12.3 ±0.5 


if* -+ ifvr 




8.6 


7.8 


7.8 ±0.2 


7.2 


7.5 ± 0.2 


7.7 ± 0.2 


7.5 ± 0.2 


7.1 ± 0.4 


4> -> if if 




2.8 


2.5 


2.5 ±0.1 


2.3 


2.3 ± 0.1 


2.4 ± 0.1 


2.3 ± 0.1 


1.9 ±0.2 


/ 2 (1275) - 


7T7T 


9.5 


8.6 


8.7 ±0.3 


9.2 


8.4 ± 0.3 


9.2 ± 0.3 


9.4 ± 0.2 


12.5 ± 1.3 


/ 2 (1275) -+ 


if if 


2.5 


2.3 


1.7 ±0.1 


2.3 


1.4 ± 0.1 


1.6 ± 0.1 


1.4 ± 0.0 


2.9 ±0.4 


/ 2 (1275) - 


m 


0.7 


0.8 


0.4 ±0.4 


0.7 


0.3 ± 0.0 


0.4 ± 0.0 


0.3 ± 0.0 


0.9 ±0.2 


/o(1300) -» 


1T7V 


11.5 


10.5 


15.5 ± 1.7 


5.5 


19.6 ± 1.6 


17.2 ± 1.6 


16.5 ± 0.4 


16.0 ±6.7 


/o(1300) - 


KK 


9.5 


8.6 


10.3 ±0.4 


6.3 


9.6 ± 0.4 


9.7 ± 0.3 


8.7 ± 0.2 


4.5 ± 2.1 


a 2 (1320) -> 


rjn 


2.0 


1.8 


1.7 ±0.1 


2.0 


1.5 ± 0.1 


1.4 ± 0.1 


1.6 ± 0.0 


3.9 ±0.4 


a 2 (1230) -> 


KK 


4.2 


3.8 


3.0 ±0.2 


3.9 


2.6 ± 0.2 


2.6 ± 0.1 


2.7 ± 0.1 


2.3 ±0.3 


/ 2 (1525) -» 


KK 


9.4 


8.5 


8.4 ±0.3 


8.1 


7.5 ± 0.2 


8.3 ± 0.3 


8.1 ± 0.2 


7.3 ± 1.6 


/ 2 (1525) -+ 




2.6 


2.4 


2.3 ±0.1 


2.4 


1.9 ± 0.1 


2.1 ± 0.1 


2.1 ± 0.1 


4.6± 1.1 


If i (1275) - 


* if O) 


7.2 


6.5 


6.1 ±0.2 


5.5 


5.4 ± 0.1 


5.3 ± 0.1 


4.8 ± 0.1 


3.1 ± 0.9 


if i (1400) - 


♦ Ku> 


1.9 


1.7 


1.9 ±0.1 


1.4 


1.7 ± 0.0 


1.6 ± 0.0 


1.5 ± 0.0 


1.3 ±0.8 


if* (1410) -» ifvr 


1.3 


1.2 


1.1 ±0.2 


0.1 


1.3 ± 0.1 


2.6 ± 0.3 


1.9 ± 0.0 


3.9 ±0.8 


if (1430) - 


* Kn 


7.5 


6.8 


10.4 ± 1.1 


3.3 


12.5 ± 1.0 


11.1 ± 1.0 


10.3 ± 0.3 


16.3 ±2.0 


if 2 (1425) - 


♦ Kir 


7.8 


7.1 


7.1 ±0.2 


7.3 


6.8 ± 0.2 


6.7 ± 0.2 


7.4 ± 0.2 


7.2 ± 1.5 


if 2 (1425) - 


♦ Ku> 


2.1 


1.9 


1.3 ±0.1 


1.8 


1.1 ± 0.1 


1.2 ± 0.1 


1.1 ± 0.0 


1.7 ±0.6 


if 2 (1425) - 


* Krj 


4.2 


3.6 


3.1 ±0.2 


4.0 


2.7 ± 0.2 


2.5 ± 0.2 


2.7 ± 0.1 


0.4 ±0.5 


if i (1680) - 


♦ Kn 


3.5 


3.2 


5.1 ±0.4 


1.5 


4.9 ± 0.4 


5.0 ± 0.3 


3.9 ± 0.1 


11.2 ± 3.0 


if 3 (1780) - 


♦ Kn 


4.8 


4.3 


4.3 ±0.2 


5.1 


4.6 ± 0.2 


3.7 ± 0.2 


5.3 ± 0.1 


5.6 ±0.8 


if 3 (1780) - 


* Krj 


3.0 


2.7 


2.1 ±0.2 


3.3 


2.0 ± 0.2 


1.4 ± 0.2 


2.2 ± 0.1 


3.6 ±0.8 



6 



enhancement of large decay amplitudes, while keeping the others at approximately constant 
values, thus increasing the agreement with the experimental data. The most noticeably 
different fit is Fit IV. 

The effective of modifying the wave function is less obvious, although the wave functions 
expanded on a larger basis (N = 3) gives slightly better fits than those expanded on the 
smaller basis. 

We have also fitted using only the first three decays of table [TIL namely p — > tttt, K* — > Kit 
and <p — > KK. We have also omitted various ones of the decays in table [H] from the fit, 
such as the decays of the / (1300). Except for a very few isolated amplitudes, none of these 
modifications have had a significant effect on the decay amplitudes that result. We note 
also that non-relativistic phase space consistently gives by far the worst results, while the 
weak-binding prescription gives the best results. However, we find that in the case of the 
full relativistic phase space, there are two decays that consistently contribute the most to 
the 'badness of fit'. One is the decay p — > 7r7r, while the other is -£f 2 (1425) — > Krj. But for 
these two decays, relativistic phase space and the weak-binding prescription would give fits 
with essentially the same values of x 2 - 

The wave functions that we use take into account SU(3) breaking effects. Strictly speak- 
ing we should also fully include such effects in the amplitudes. These effects arise from three 
sources. The most obvious of these is in the calculation of phase space, where we use either 
the physical masses of the states (for non-relativistic and relativistic phase space), or the 
'weak-binding' masses appropriate to the mesons in the decay. The second source of SU(3) 
breaking is due to the fact that the amplitude depends on the momentum of the daughter 
mesons, calculated in the rest frame of the parent. Here we also use the physical masses 
of the states to calculate this momentum. The third source of SU(3) breaking arises in the 
evaluation of the matrix element of the transition operator. The masses of the quarks or, 
more precisely, various ratios of masses of the quarks, enter explicitly into this calculation. 
We have found that treating the u, d and s as degenerate provides the best results. Thus, 
SU(3) breaking at this level is not necessary. 

A few comments on the results of table [TJ are in order. The amplitudes for the i^i(1275) 
and i^i(1400) are calculated from those for the 3 Pi and 1 P\ states of the model spectrum, 
assuming a mixing angle of 34°. In tables |T| to [V], we ignore this mixing and present results 



for the unmixed 3 Pi and l P\. The r? and rj' are assumed to be 

1 



1 

V 



V2 



(uu + dd) — ss 



V2 



" = 7i 



-j= (uu + dd 



+ ss 



(11) 



respectively. The u is assumed to be purely uu + dd, and the <fi is taken as pure ss. Finally, 
in all of the results we have ignored mixings of the type 3 Lj ^ 3 L'j, while in tables |TTT| to 
0, we have also ignored mixings of the type 1 Lj <-> 3 Lj. 

The numbers in table |H] show that the model reproduces the experimental data reason- 
ably well, for all of the scenarios presented. Not surprisingly, the three-parameter fits are 
somewhat better than any of the one-parameter fits. With few exceptions, amplitudes that 
are found to be large experimentally, also turn out to be large in the model, and the same 
is generally true of small amplitudes. The few exceptions where a large discrepancy exists 
are /o — > KK, X 1 (1275) — > Kuj and i^i(1680) — > Kir. In the case of the first of these 
decays, we note simply that the status of this state is very much in question, so that the 
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TABLE III. Results from various fits, for states with 1 = 1. The key is as in table 0. 
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a Experimental values are for the sum over all relevant partial waves. 

partial widths themselves are also uncertain. In the case of the i£i(1275), we note that the 
mixing we have used has not resulted from a 'rigorous' calculation, but is simply 'borrowed' 
from the literature. In the case of the i^i(1680), there is the possibility of mixing with the 
radially excited 3 Si state, not taken into account in this calculation. However, indications 
are that this mixing is very small 0, and would probably not play a large role in bringing 
the model calculation into closer agreement with the experimental measurement. 

In tables [TTT| to M, we have extended the model calculation to examine decays of a few, 
semi-randomly selected, radially and orbitally excited states. For these tables, the results 
correspond to the first three fits of flj, respectively. We preface our discussion here by 
reminding the reader that the assignment of many of the experimentally observed states is 
still not clear, and it is not our aim, at least not in this work, to clarify the assignments of 
any of these states, but simply to see how well our model for their decays (and implicitly, of 
their spectrum) works. Thus, we have chosen only a few of the 'better established' states for 
this discussion. Note that if the predicted value for an amplitude is less than 0.05 MeV 1 / 2 , 
we do not show the amplitude. 

In obtaining the numbers for decays with a broad vector meson (p or K*\ the uj is less 
than 10 MeV wide, and we treat it as a narrow state in our calculation) in the final state, 
we integrate over the line shape of the broad final state. That is, for such decays, the decay 
rate is taken as 

iw = r~ dk W)l a rv(fc) 

{mA -E B (k)-E v (k)) 2 + T -^' 
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TABLE IV. Results from various fits, for states with 1 = 0. The key is as in table [n|. 





decay 




pl/2 


(MeV 1 / 2 ) 








Fit I 


Fit II 


Fit III 


experiment 


/a (1275) -» 


7T7T 


9.5± 0.1 


8.6± 0.1 


8.7± 0.2 


12.5 ± 1.3 




KK 


2.5± 0.1 


2.3± 0.1 


1.7± 0.1 


2.9 ±0.4 




r\r\ 


0.7± 0.1 


0.8± 0.1 


0.4± 0.1 


0.9 ±0.2 


fi(1280) — ► 


KK* 


0.5± 0.1 


0.4± 0.1 


0.4± 0.1 


not seen 


fii (13001 — » 


7T7T 


11. 6± 5.7 


10. 5± 5.1 


15. 6± 6.6 


16.0 ± 6.7 




/s: A' 


9.5± 3.6 


8.6± 3.2 


10.3± 2.5 


4.5 ± 2.1 




rjrj 


6.1± 6.1 


5.5± 5.5 


6.2± 6.2 


seen 


fof 15251 — > 


KK 


9.4± 0.2 


8.5± 0.2 


8.4± 0.2 


7.3 ± 1.6 




an 


2.6± 0.1 


2.4± 0.1 


2.3± 0.1 


4.6 ± 1.1 




vv' 


0.0± 0.1 


0.0± 0.1 


0.0± 0.0 






KK* 


3.2± 0.5 


2.9± 0.5 


2.2± 0.4 




/i(1510) -» 


(KK*) L=0 


9.4± 0.5 


8.5± 0.4 


8.5± 0.5 


seen a 




(KK*) L=2 


1.9± 0.3 


1.7± 0.3 


1.3± 0.2 




/o(1500) -> 


KK 


12. 6± 0.3 


11. 4± 0.3 


15.5± 0.3 






VV 


7.3± 0.1 


6.6± 0.1 


8.5± 0.1 


seen 




rm' 


0.0± 3.9 


0.0± 3.6 


0.0± 3.6 


seen 


0(1680) -> 


KK 


5.6± 0.2 


5.0± 0.2 


6.8± 0.1 


seen 




m 


1.5± 0.1 


1.4± 0.1 


1.8± 0.0 






rm' 


1.4± 1.4 


1.2± 1.2 


0.8± 0.8 






KK* 


12. 6± 1.9 


11. 4± 1.7 


11. 2± 2.0 


dominant 



Experimental values are for the sum over all relevant partial waves. 

where V denotes the vector meson, and Ty is the total width of the vector meson. The 
errors on the theoretical numbers in these tables arise from evaluating the amplitudes at the 
upper and lower limits of the mass ranges allowed by the experimental error on the masses 
of the states. 

As with the fitted decays, we find that the overall agreement of the model predictions 
with the experimental data is quite good. There are a few instances in which there is 
noticeable disagreement, and we comment on these below. 

Our results for the ai(1235) are 'consistent' with experiment, in that the decay of this 
state to pir is 'dominant' (of a total width of about 400 MeV), while we predict a pir partial 
width of about 100 MeV. Not many other exclusive final states have been identified in the 
decays of this meson, leaving the possibility that much of its width comes from decays to 
excited mesons (or exotics). For the rest of the decays of the 1=1 states, the results that 
we obtain are in reasonable agreement with what little data there is. 

In the 1 = channels, our predictions are in quite good agreement with experiment, 
except for the decays of the /b(1300). However, the assignment of this state is so uncertain, 
(e. g.; there still is the question of whether this is a single state or not) that interpretation of 
the discrepancy is difficult. In this sector there is one decay that is worthy of some comment. 
At its nominal mass, the /o(l°00) will not decay to rji]', yet the partial width for this decay 
has been listed as 'large (dominant)' JT3]. However, we see that if we change the mass of 
this state by 10 MeV in our model (a small amount, considering how broad ++ states are 
expected to be), the amplitude for decay to rjrj' grows rapidly. This is not surprising, as 
this is an S"-wave amplitude, which grows roug hly as (M - M th y/ 2 , where M is the mass of 
parent and M th is the threshold for the decay. This would suggest that either (a) this state 



is more massive than reported in flI2|; or (b) the amplitude grows so rapidly that we should 
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integrate over the line-shape of at least the rf in calculating the amplitude. This decay also 
illustrates some of the importance of threshold effects. 

In the kaon sector, our predictions are again in reasonable agreement with experiment. 
In the case of the KiS, we can not really compare theory and experiment, since we have not 
included mixing. This can easily be taken into account in the decays to the narrow states, 
but not for the broad states, since an integral over the partial width is required to evaluate 
the amplitude. Thus, simply constructing a linear combination of the reported amplitudes 
will be ignoring potentially large interference effects. 

The place where our model 'fails' most noticeably is in describing the decays of the 
K 1 (1680). There, our amplitudes are consistently smaller than the experimental ones. Note, 
however, that for the three measured amplitudes, our results agree 'qualitatively' with the 
experiments, in that the 3 amplitudes are all predicted to be of roughly the same size. 

One other aspect of our results requires some comment. Our results have been largely 
independent of the size of the expansion basis used, as well as of the choice of potential 
used to generate the wave functions. This may be understood by examining what part of 
the wave function provides the dominant contribution to any one of the amplitudes. It is 
apparent that these decays are neither fully 'long-distance' nor 'short-distance' processes, 
but more like 'intermediate-distance' (1-5 GeV -1 ) processes. Changes in the basis size have 
their main effect on the shape of the wave function in the region of the origin, as well as 
on the long distance tail. By changing the potential from All to API, we change the same 
regions of the wave function. 

IV. CONCLUSION 

We have examined the predictions of a quark model for the strong decays of several 
mesons. While this study was not exhaustive, it has been enough to show that the model 
enjoys some success, but that some improvements may be possible. The main improvements 
that can be made exist at the level of the meson wave functions. 

At present, there is no spin-orbit term in either of the potentials used, nor is there a 
tensor interaction. Both of these will lead to modifications of the wave functions that will, 
in turn, lead to changes in the strong decay amplitudes predicted by the model. These 
terms will also induce various mixings among the states, which will also modify the wave 
functions. Such a study is underway, and the use of the modified wave functions in a strong 
decay analysis is left as a possible future endeavor. It would certainly be of interest to 
apply such wave functions to a more exhaustive examination of the strong decays of mesons, 
as there are several puzzles outstanding in the meson sector, particularly for mesons with 
masses between 1.0 and 2.0 GeV. To date, models of this kind have been our main source of 
insight into the meson spectrum. Barring a radical revolution in the near future, they will 
remain so, at least for some time to come. 
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TABLE V. Results from various fits, for kaons. The key is as in table [TI| 





decay 




F 1 /' 2 (MeV 1 ^) 


| 






Fit I 


Fit II 


Fit III 


experiment 


Ki( L Pi) - 


Kw 


0.0± 2.9 


n n_l_ ft 


n n_l_ 


3.1 ±0.9 a b 




(Kp) L=0 


3.0± 0.3 


2.7± 0.3 


2.5± 0.3 


6.1 ± 1.1 a b 




{Kp) L=2 


0.9± 0.1 


0.8± 0.1 


0.6± 0.1 






(irK*) L=0 


4.7± 0.1 


4.2± 0.1 


4.2± 0.1 


3.8 ±1.0"' 




(nK*) L=2 


1.8± 0.1 


1 ft.4- n 1 


1 ix n 1 

± . 1 _i_ U.l 




#i( 3 Pi) -* 


(Kw) L=0 


6.9± 0.1 


6.2± 0.1 


6.1± 0.0 


1.3 ±0.8"' 




(Kw) L=2 


1.3± 0.1 


1 94- n 1 
l.zzt U.l 


n 74- n 1 
U. 1 ± U.l 






(K P ) L =o 


8.5± 0.2 


7.7± 0.2 


7.7± 0.3 


2.3 ±1.2°' 




(Kp) L =2 


2.4± 0.2 


2.2± 0.2 


1.6± 0.1 






(itK*) l=0 


10.3± 0.1 


9.3± 0.1 


10. 5± 0.2 


12.8 ± 0.9 a b 




(nK*) L=2 


4.5± 0.2 


4.1± 0.2 


3.1± 0.2 




K*(1410) — > 


Kit 


1.3± 0.1 


1.2± 0.1 


l.lrt 0.0 


3.9 ± 0.8 6 




Kr) 


0.1± 0.0 


0.1± 0.0 


0.8± 0.0 






Kw 


2.8± 0.1 


2.5± 0.1 


1.5± 0.1 






Ko 


3.8± 0.1 


3.4± 0.1 


2.3± 0.1 


< 4.0 b 




ttK* 


4.9± 0.1 


4.4± 0.0 


3.3± 0.1 


> 9.5 6 


ATo(1430) -» 


Kit 


7.5± 0.2 


6.8± 0.2 


10. 4± 0.3 


16.3 ± 2.0 




Kn 


11. 2± 0.1 


10. 2± 0.1 


13. 2± 0.1 




Xq (14251 — > 

jj V / 


Kn 


7.9± 0.1 


7.1± 0.1 


7.1± 0.2 


7.2 ± 1.5 




Kn 


4.2± 0.1 


3.7± 0.1 


3.1± 0.1 


0.4 ± 0.4 




Kw 


2.1± 0.2 


1.9± 0.2 


1.3± 0.1 


1.7 ± 0.6 




Kn 

1 


3.7± 0.2 


3.3± 0.2 


2.5± 0.2 


3.1 ± 0.2 




wK* 


6.4± 0.1 


5.8± 0.1 


4.4± 0.1 


5.2 ± 0.3 


K(1460) -> 


Kw 


5.8± 0.7 


5.2± 0.6 


4.1± 0.6 






Ko 


8.7± 0.8 


7.8± 0.7 


6.6± 0.8 


~ 5.8 




■kK* 


12. 2± 0.9 


11. 0± 0.8 


10. 1± 0.5 


~ 10.1 


ffoflSROl — » 


(Kw)t 1 


2.4± 0.2 


2.2± 0.2 


2.1± 0.5 






(Kw) t ■» 


1.5± 0.9 


1.3± 0.8 


0.6± 0.7 






(Kn\r 1 


4.4± 1.5 


4.0± 1.3 


4.0± 1.6 






(Ko)t-% 


3.1± 1.7 


2.8± 1.5 


1.6± 1.0 






(nK*\r 1 


3.0± 0.5 


2.7± 0.4 


2.8± 0.9 


seen b 




(■kK*)l=3 


2.3± 2.0 


2.1± 1.8 


1.2± 1.5 






(ljK* ) Q— 1 r — 1 


0.0± 0.8 


0.0± 0.7 


0.0± 0.4 






(u)K* ) Q 1 T — "i 


0.0± 0.1 


0.0± 0.1 


0.0± 0.0 




A'i(1680) — ► 


Kn 


3.6± 0.1 


3.2± 0.1 


5.2± 0.1 


11.2 ± 3.0 




Kn 


4.0± 0.1 


3.7± 0.0 


5.0± 0.0 






Kuj 


2.0± 0.1 


1.8± 0.1 


2.1± 0.0 






Ko 


3.4± 0.1 


3.1± 0.1 


3.6± 0.0 


10.1 ± 2.5 




TTK* 


3.2± 0.2 


2.9± 0.1 


3.8± 0.3 


9.8 ± 0.9 




(ujK*)Q—r\ T— 1 


1.2± 0.6 


1.1± 0.5 


0.7± 0.3 






(ujK* ) q 1 /" 1 


0.5± 0.3 


0.5± 0.2 


0.3± 0.1 






((jjK* ) Q 7" ^ 


0.4± 0.2 


0.4± 0.2 


0.1± 0.1 




Ko(l 77(1) — > 


(Kt,]\r 1 
UJ )]_, — 1 


2.3± 0.1 


2.1± 0.1 


2.8± 0.0 


seen b 






2.6± 0.1 


2.4± 0.1 


1.7± 0.1 








4 3± 0.0 


3.9± 0.0 


4.9± 0.0 






(Ko\t Q 


5.1± 0.3 


4.6± 0.3 


3.6± 0.3 






(ttK*)t i 


4.9± 0.1 


4.5± 0.1 


6.3± 0.2 






(vrA*) L =3 


6.8± 0.5 


6.2± 0.4 


5.0± 0.4 






(ojA"*)s = 2,l=1 


3.5± 0.7 


3.2± 0.6 


2.2± 0.4 






(a;A"*) s=2 ,L=3 


1.5± 0.4 


1.4± 0.3 


0.5± 0.1 




A 3 (1770) -+ 


Xtt 


4.8± 0.1 


4.3± 0.1 


4.3± 0.1 


5.6 ±0.8 




Kn 


2.9± 0.1 


2.7± 0.0 


2.1± 0.1 


3.6 ± 0.8 




Kw 


2.9± 0.1 


2.6± 0.1 


1.9± 0.1 


3.6 ± 0.5 




Kp 


6.0± 0.3 


5.4± 0.3 


4.3± 0.3 


8.6 ± 0.8 




nK* 


8.0± 0.5 


7.2± 0.5 


5.8± 0.5 


6.7 ± 0.7 




(wK*) s=0:L=3 


0.4± 0.1 


0.4± 0.1 


0.2± 0.0 






{wK*)s=2,L=l 


6.7± 1.5 


6.1± 1.4 


4.3± 1.0 






(wK*)s=2,L=3 


1.0± 0.3 


0.9± 0.2 


0.3± 0.1 





a The experimental numbers are for the fully mixed states, while the theoretical numbers for the 
1 P\ and 3 Pi unmixed states. 

Experimental values are for the sum over all relevant partial waves. 
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